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ABSTRACT 



The single-blow transient testing facility of the Naval Postgraduate 
School was reconstructed. With this apparatus the flow friction and 
heat transfer characteristics of four grades of steel wool, considered 
to be representative of fibrous metals, were determined. The flow 
friction characteristics were determined on the basis of the isothermal 
pressure drop across the medium, while the heat transfer results were 
determined by means of the single-blow transient testing technique. 

The heat transfer results are presented as the Colburn j-factor. 

The flow friction results are presented in two forms: as the Fanning 

friction factor, f; and as the correlating friction factor, C^. All 
results are presented as a function of Reynolds number based on the 
product of the experimentally determined permeability and tortuosity 
parameter of each sample. 
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NOMENCLATURE 



English Letter S 3 mibols 



Typical Units 



a Linear coefficient of Reynold's 

equation, reciprocal of the 
permeability (1/K) 

A Matrix total heat transfer area 



ft 



-2 



ft 



2 



A Matrix minimum free flow area 

c 

Aj. Matrix total frontal area 

fr 

A Matrix solid cross-sectional area 

g 

available for longitudinal conduction 



ft 



2 



ft 



2 



ft 



2 



b Quadratic coefficient of Reynold's ft 

equation, interpreted as a tortuosity 
parameter 



c Cubic coefficient of the empirical 
equation for the pressure gradient 



dimensionless 




c 

P 

c 

s 

C 

s 

d 



Correlating friction factor 

Fluid specific heat at constant pressure 

Matrix material specific heat 

Matrix thermal capacity (W c ) 

s s 

Rod diameter in stacked screen model 



dimensionless 

BTU/(lbm-°F) 

BTU/(lbm-°F) 

BTU/°F 

ft 



djj Fiber hydraulic diameter 

D Flow passage hydraulic diameter 

E Friction power per unit area 

G Mass flux (m/A ) or (m/A ) 

c fr 

g Gravitational acceleration 

g Proportionality factor in Newton's 

Second Law 



ft 

ft 

hp/ft^ 

lbm/(hr-ft^) 

ft/sec^ 

32.2 lbm-ft/(lbf-sec^) 



2 

h Surface heat transfer coefficient for BTU/(hr-ft 

convection; heat transfer power per 
unit area per degree temperature dif- 
ference 
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Fluid thermal conductivity 


BTU/(hr-ft-°F) 


Matrix permeability 


ft^ 


Contraction coefficient 


dimensionless 


Expansion coefficient 


dimensionless 


Matrix thermal conductivity 


BTU/(hr-ft-®F) 


Length parameter in 


ft 


Total matrix flow length 


ft 


Mass flow rate 


Ibm/hr 


Pressure 


Ibf/ft^ 


Matrix porosity (l-W^/(yO^A^^L) ) 


dimensionless 


Heat transfer rate 


BTU/hr 


Specific gas constant (53.35 for air) 


ft-lbf/(lbm-®R) 


Flow passage hydraulic radius (D /4) 

n 


ft 


Temperature 


OF 


Absolute temperature 


OR 


Normalized temperature 


dimensionless 


Flow velocity 


ft/sec 


Bulk velocity (G^) 


ft/sec 


Specific volume (1/^) 


ft^/lbm 


Matrix volume 


ft^ 


Average fiber width 


in 


Fluid mass in matrix 


Ibm 


Matrix mass 


Ibm 


Uncertainty interval for quantity y 


units of y 


Distance along flow passage from the 
matrix inlet 


ft 


Pitch of rods in stacked screens model 


dimensionless 


Reduced length (N^^(x/L)) 


dimensionless 



Greek Letter Symbols 



dimensionless 

sec, hr 

lbm/(hr~ft ) 
Ibf-hr/ft2 

dimensionless 

Ibm/ft^ 

ft^/hr 

Siibscripts 

atm Local atmosphere 

ave Average 

f Fluid (gas, air) 

i Initial; inlet 

m Matrix; mean 

0 At orifice 

s Solid (matrix material) 

STD Standard (temperature and pressure) 

X Local conditions 

1 Inlet conditions (upstream of matrix and heaters) 

2 Inlet conditions at matrix entrance 

3 Exit conditions at matrix outlet 

Dimensionless Groupings 

Correlating friction factor; ratio of dissipative 
effects to inertial effects 



/« 

A 

e 

t 



Specific surface ( compactness A/A^^/L ) 

Ratio of orifice diameter to pipe 
diameter 

Difference or change (of time, tern- 
perature, etc.) 

Time 

Fluid viscosity 



Free time 
Density 

Kinematic viscosity 
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f 

j 

A 



Fanning friction factor; ratio of wall shear stress 
to fluid dynamic head 



2/3 ■ 

Colburn j factor )• This factor plotted versus 

Reynolds number defines the surface heat transfer char- 
acteristics 



Longitudinal heat conduction parameter for solid material 

(k A /mLc^ . ,) 
s s fluid 



‘ Time parameter (hA0/W^c^) 

Np^ Prandtl number; ratio of viscous dissipation to thermal 

dissipation (i.e., kinematic viscosity /thermal diffusivity) 

N Reynolds number; ratio of inertial effects to viscous 
® effects (^Ul^) 

Ng^ Stanton number; ratio of heat transfer power per unit 
surface area to fluid thermal capacity per unit flow 

area (h/Gc ) 

P 

Number of transfer units; ratio of total heat transfer 

TU 

power to total fluid thermal capacity (hA/mc ). 
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I. INTRODUCTION 



Regenerators add on the order of ten percent to the thermal effi- 
ciency of the basic gas turbine power cycle. Due to the current 
emphasis on the use of gas turbines, as well as the reawakening interest 
in the Stirling and Ericsson cycles, where a regenerator is an integral 
part of the cycle, much attention has been paid to the design of these 
heat exchangers. Frequently the applications of such cycles are weight 
or volume-limited, and therefore regenerators must possess a high heat 
transfer surface density. 

However, as noted by Kays and London in their monograph on compact 
heat exchangers [13]^, increasing the heat transfer area necessarily 
increases the pressure drop (hence the pump work required) across the 
heat exchanger, particularly when the working fluid is a gas. Thus, 
the consideration of the heat transfer characteristics of a surface for 
compact heat exchanger applications cannot be divorced from a consid- 
eration of the friction characteristics of the surface. 

Between 1965 and 1968 several compact heat exchanger surfaces were 
extensively studied at the Naval Postgraduate School. These included 
perforated nickel plate matrices and offset rectangular fin matrices 
of the AiResearch Manufacturing Company studied ^by Piersall [23]; six 
different plate-and-fin cores by the Solar Division of International 
Harvester studied by Ball [2], Traister [29], and Trost [30]; and the 
ceramic cores by Corning Glass, studied by Bruno [5]. All of these 



Numbers in brackets refer to items in the bibliography. 
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surfaces were fundamentally porous to flow in one direction only and 
were intended for use in gas turbine plants. 

Fibrous metals (such as steel wool) possess the high surface density 
requisite for a compact heat exchanger. Moreover, the apparent small 
hydraulic diameter of the effective flow channels of such materials 
should lead to large values of the convective heat transfer coefficient 
h, since h always varies with a negative power of the channel hydraulic 
diameter [13]. Also the tortuous nature of the three-dimensional flow 
through fibrous media hinders the development of the momentum and thermal 
boundary layers, which should further increase the value of h. Although 
the high surface area to volume ratio of the individual fibers excludes 
fibrous metals from high temperature environments (e.g. , the open cycle 
gas turbine regenerator), such materials do have application in low- 
temperature closed power cycles and in refrigeration cycles such as the 
reversed Brayton, reversed Stirling, Gifford-McMahon , or Pulse Tube 
cycles . 

This thesis was directed at four objectives: 

1. Reconstruct the transient test apparatus (dismantled since 
1968), including new instrumentation. 

2. Test a previously studied compact heat exchanger surface 
to insure reliability of the reconstructed apparatus. 

3. Determine the heat transfer and flow friction character- 
istics of various grades of fibrous metals, specifically 
steel wool. (The ”grade” of a fibrous metal refers to 
the size of the fiber, and not necessarily to the porosity 
of the fibrous metal matrix. Porosity was not an inde- 
pendent test variable). 

4. Ascertain the appropriate length parameter to be used in 
the definition of a Reynolds number upon which a unified 
correlation of the heat transfer and flow friction char- 
acteristics may be based. 
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II. SUMMARY OF -THEORY 



A. REYNOLDS NUMBERS FOR POROUS MEDIA FLOWS 

Flovr friction and heat transfer data are usually correlated with 
the flow Reynolds number for the surface geometry under consideration: 



where 1 is a characteristic length parameter for the surface geometry. 
The results thus correlated are independent of any particular geomet- 
rically similar surface. For flow in fibrous media the dilemma is the 
selection of the characteristic length parameter. If the fibers were 
packed so tightly that the flow could be unambiguously considered 
internal in nature , the obvious choice would be some sort of effective 
hydraulic diameter of a representative flow channel. This was the 
procedure in the previously tested cores (e.g. , the Solar series, 
Cercor, etc.), which were porous to flow in one dimension only and 
possessed a well-defined channel hydraulic diameter. At the other 
extreme, if only one fiber of the medium were under consideration 
(i.e., purely external flow), the characteristic length parameter 
would be the hydraulic diameter of the fiber. However, neither such 
simple case exists. 

Three approaches were made to the selection of the characteristic 
length parameter. In the first, fibrous metals were assumed to be 
porous to flow in only one dimension and the equations applied to the 
previously tested cores were used. The second approach is based on 
the in-line stacking of crossed-rod matrices. Both of these idealized 
models produce a physically conceivable characteristic length. The 

third approach was based on the fundamentals of flow in porous media, 
but produces no such conceivable dimension. 
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1. Uni-dimensional Porosity Model 



Consider a heat exchanger with n uniform channels of length L. 

If A is the free flow area of the entire cross-section, and A is the 

heat transfer area of the entire exchanger, then: 

Free flow area/channel = A /n 

c 

Nominal volume /channel = (A /n)L 

c 

Interior surface/channel = A/n 

But the interior surface of a channel equals the cross-sectional perim- 
eter times the length of the channel. Therefore 
Perimeter /channel = A/(nL) 

The hydraulic diameter of the channel is defined as four times the cross- 
sectional area of the channel divided by the wetted perimeter, or: 



X) = 4 — ^ = 4 

^ (a/[ol]^ a 



( 1 ) 



For such a model which is porous to flow in only one dimension, the 
porosity (i.e., volume of voids over nominal volume) is: 






where equals the total frontal area of the medium. Since the 
nominal volume of the exchanger is the product of and L, the 
specific surface (i.e. , heat transfer area per unit of volume) is, 
with substitution from (1) and (2): 



A = z 4 ^ = 4 ^ 



(3) 



\ 
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As explained in Appendix B, the porosity and heat transfer area (hence 
also the specific surface) may be determined independently of these 
relations for fibrous metals. Thus the ‘equivalent channel hydraulic 
diameter for this model may be computed from (.3). This is the charac- 
teristic length used in the Reynolds number under this idealization. 

2. In-line Crossed Rod Matrices Model 

This model consists of square-mesh cross-rod screens which are 
stacked so that there is no separation between layers and so that the 
mesh in one layer coincides with the meshes in adjacent layers (see 
inset). The rods are assumed to be of diameter d and separated by 
uniform distance x^d. 

On the basis of a repeatable unit cell (similar to the procedure 
in the preceding section) the following relations are easily derived [13] 

^ = 1 - 1^/(4 (4) 

Ac/Af,=^ [K- <4a) 

^ = Tr/(Ktd^ (5) 

Dh = 4k^/tt - 1 (6) 

The pitch, x^, is calculated from the known porosity, with the use 
of (4). Again this equivalent channel hydraulic diameter is used as the 
length parameter for this model. Note that this model differs from the 
previous in that it is porous to flow in three dimensions, rather than 
one. For the cases where the fiber is not of circular cross-section, 
the rod diameter d is replaced by the fiber hydraulic diameter in the 
equation for the specific surface (5). 







K ,<1 
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3. Flow in Porous Media 

Analysis of the fundamentals of flow in porous media provides 
a length parameter which is more justifiable, in the sense that no 
geometrical approximations are made. 



sand filters of the water purification plant at Dijon, France, made the 
empirical observation that the bulk (i.e., nominal) flow rate of fluid 
through the filters is proportional to the difference between the fluid 
inlet and outlet heads, and inversely proportional to the thickness of 
the filter. Since this observation was valid only macroscopically , 
and therefore over finite distances, it was later formulated on a dif- 
ferential basis as what is now known as Darcy’s Law of fluid flow in 
porous media: 



Here U is the bulk velocity of the fluid flow (i.e. , based on the total 
frontal area of the medium and not the internal free flow area), p is 
the fluid pressure in the medium, and g is the gravitational accelera- 
tion. The constant K is known as the permeability of the medium. 



of fluids in porous media. The significance of this law is that it 
relates the internal pressure field to the bulk velocity U, a quantity 
which is measured external to the medium. Together with the continuity 
equation for fluid flow in a porous mediiim: 



In 1856 Darcy, after a study of the flow of water through the 




(7) 



Darcy’s Law is the fundamental relation describing the flow 




( 8 ) 



total porosity of the medium). 
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Darcy’s Law completely specifies the flow field. Under steady state 
conditions these two equations then reduce to Laplace’s equation, and 
a solution is possible. It must be emphasized that Darcy’s Law is an 
empirical relation; it has no analytical basis. Theoretically, if a 
satisfactory model can be formulated for the porous medium, and suf- 
ficient boundary conditions stated, the same result could be achieved 
by a. solution (closed form or numerical) of the Navier-Stokes equations 
(which would involve the internal velocities rather than the external, 
bulk velocity). In practice this analytical approach has rarely been 
successful, largely because of difficulties with the boundary conditions. 

It is convenient at this point to define some terms of impor- 
tance in the study of Darcy’s Law: 

a) Porosity,^, is the volume of voids divided by the 
total (i.e., nominal) volume. 

b) Permeability, K is , a measure of the ease with which 
a fluid may transverse a medium under the influence 
of a given pressure gradient. It is important to 
note that the permeability has units of length squared 
and is independent of the fluid. Permeability is 
defined by Darcy’s Law (7). 

c) Tortuosity is the relative average length of the 
flow path of a fluid particle from one side of the 
porous medium to the other. 

Clearly the permeability is the important parameter of the 
medium for predicting the pressure drop across it. Unfortunately, it 
is the porosity of the medium, and not the permeability, which is 
generally known. For a given porous medium, the permeability is a 
unique, but unknown function of the porosity [22, 27]. Several theories 
have been advanced which purport to theoretically relate porosity and 
permeability. The oldest is that due to Kozeny (1927) but his own 
experimental data shows discrepancies from the theory of -68.7 to +85.7 
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per cent [22]. The most recent theory is based on statistical mechanics, 
but regretably attaches the randomness to the fluid flow, and not the 
medium, where it properly belongs [27]. Thus it may be seen that the 
most accurate determination of the permeability is by experiment. 

If Darcy’s Law C7) is written for one-dimensional flow in the 
absence of significant changes in gravitational potential, it becomes 
(after rearrangement): 

dp _ (9) 

dx ■ K 

Obviously the permeability, K, can be determined by measuring the iso- 
thermal pressure drop across a thin slice of medium at different bulk 
velocities, U. (Constant temperature is required so that jx remains 
constant) . 

Mother Nature has not been so kind, however, for Darcy’s Law is 
strictly valid only for low flow rates. At higher velocities the 
magnitude of the pressure gradient increases with velocity more rapidly 
than the linear rate indicated by Darcy’s Law. Based on experimental 
data gathered in the non-Darcy regime, Reynolds (1900) suggested the 
relation: 

o.iuU + (10) 

F 3c 

Forchheimer (1901) added a cubic term to better fit the data. However, 
the coefficient of the cubic term was found to be of such small magnitude 
that the term is generally neglected. Missbach (1937) proposed that the 
pressure gradient in non-Darcy flow be proportional to U to some exponent 
between one and two. As this relation could not even be partially jus- 
tified theoretically, it has generally been ignored in favor of Reynold’s 
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equation. The most refined formulation is that of Polubarinova-Kochina, 
which includes momentum changes due to compressibility effects : 




( 11 ) 



This additional term is the acceleration across the medium due to the 
decrease in the density of the fluid. Since both a stationary frame of 
reference and a frame moving at the bulk velocity are Inertial frames, 
this time-dependent acceleration may be rewritten as a convective accel- 
eration [8]: 



tified on quasi-theoretical grounds by analogy to the pressure drop in 
pipe flow. In laminar pipe flow the pressure gradient is proportional 
to the velocity (the Hagen-Poiseuille equation), while in turbulent pipe 
flow the gradient is proportional to the square of the velocity. Since 
transition between laminar and turbulent flow in fibrous media would not 
be as sharply defined as in pipe flow, the reasoning went, the pressure 



loss must be the sum of the two analagous losses [22]. 

It will be recalled that the non-linear (viz. , inertial) terms 



Reynolds number exceeds some critical value [28]. This critical value 
marks the end of the "slow flow" regime. Separation of the boundary 
layer in external flow over bluff bodies and similar effects in curved 




(11a) 




(11b) 



Although (10) is of experimental origin, it was initially jus- 



in the Navier-Stokes equations do not become significant until the - 
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channel flow occur in this new regime. For both internal and external 
flows, there is a dramatic increase in the pressure drop long before 
the flow becomes turbulent itself. It is now thought that non-Darcy 
flow is due to the emerging influence of these inertial terms [9, 26]. 
Hubbert Cl2] proposed such a completely laminar theory wherein the 
increasingly important inertial forces (which increase with velocity 
more rapidly than do the viscous forces) distort the streamlines of 
flow. This leads to increased velocity gradients and therefore, in- 
creased negative pressure gradients (by Bernoulli’s Law). These 
increased velocities also lead to increased viscous forces. The effects 
are lumped as inertial resistance, hence a greater pressure drop than 
that predicted by Darcy’s Law. 

Dimensional analysis of Reynolds equation 




2 

shows that ^ must have dimensions of one over length squared (l/L ), while 
^ has dimensions of one over length (l/L). In fact, for small bulk ve- 
locity, U, the . quadratic (i.e., inertial) term in Reynolds equation 
becomes negligible, and (10) reduces to the one-dimensional Darcy’s Law 
(9). Therefore the constant ^ must be one over the permeability (1/K). 
Thus (10) becomes: 



cjp mU , 

'dx 



( 12 ) 



The coefficient ^ of the inertial term is an inertial resistance coef- 
ficient thought to be ude to the successive contractions and expansions 
of the flow in porous media. Green and Duwez [9] assert that ^ is 
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actually a measure of the tortuosity of the medium. However, there is 
no purely theoretical reason for assuming that ^ is independent of the 
fluid [3]. 



If (10) is arranged in the following form: 



r 




l[-4p]= -L +b ^ 



(13) 



it is seen that a plot of -jj 



should be linear 



with slope equal to the tortuosity parameter, b, and intercept equal to 
one over the permeability, 1/K. Beavers and Sparrow [3] have examined 
the flow of water through a fibrous metal product known as FOAMETAL, 
which is structurally similar but not identical to steel wool. Their 
investigation covered a velocity range of 20:1 for 5 samples and showed 
no deviation in excess of 2.5% from the linear plot predicted by (13). 
This is considered to be a strong indication of the validity of Reynolds 
equation and suggests a simple means of determining the permeability and 
tortuosity parameter for a porous medium. 



taps is assumed). This approximation may be eliminated by the integra- 



under the assumption of isothermal flow of a perfect gas. Since (11b) is 
the more general of the two, it is used for the illustration of the inte- 
gration procedure. 




However, the use of this technique requires the approximation of 
(in other words, a uniform gradient between pressure 



tion of Reynolds equation (10) or the Polubarinova-Kochina equation (11b) 




(11b) 
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The viscosity is the absolute viscosity in units of force-time /length 
squared. Multiply Cllb) through by the density 







(14) 



of these identities in (14) yields: 



. _ aaQ + -tG*+ i 

9c^ 9c dxV/=/ 



r 



(15) 



The mass flux, G, is constant across the medium since it is based on the 

frontal area, A. (i.e., G = m/A^. ). Therefore 
’ fr fr 



- ^ r ^ fr {wr) 



(16) 



-R-^= ayuG 4- hG -f — p 
RTdx / 3^ 3^ 






(17) 



.-R^ = aMG -hkG’'+ p(-l/p^)4^ 



RT 



(18) 






kG'- 5!(ida 

3c 3c. ^ P dx 



_p_ dp = cy^Gdx h^dx— 



RT 



9c- V P 



(19) 



(20) 



Equation (20) can be integrated from entry to exit of the medium (a length 

L): 
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7 . 



._e 1 =ajuGL +kG^L 
2RT , r 3c 





L. 4- 



3‘ 



( 21 ) 



( 22 ) 



This equation may also be arranged so that ^the slope is again equal to 
the tortuosity parameter b while the intercept is one over the permea- 
bility (a = 1/K) : 




P -fi- - G |„( 




ZRT 3c ' 






(23) 



The choice between (13) or (23) to determine the permeability 
depends on the conditions of the experiment. The former is restricted 
to isothermal flow with a constant pressure gradient between pressure 
measurement stations. The latter is restricted to isothermal flow of 
a perfect gas, but places no further restriction on the nature of the 
pressure gradient. 

Once the permeability has been determined for a medium, a 
Reynolds number can be defined. Since is the ratio of inertial 
effects to viscous effects, and since all the viscous effects are 
embodied in the linear term in (10) or (11), while all the inertial 
effects are in the quadratic term: 



N 



Re 




/3 c 

/-WK 




.(24) 



24 



Since jx was previously identified as the absolute viscosity in force 
units 5 is the absolute viscosity in mass units. To avoid ambiguity, 

/-'r 



substitute the kinematic viscosity P for . Therefore, 



N,. = ^ (Kb) 



(25) 



It is seen that the characteristic length parameter based on this ap- 
proach is the product of the permeability and the tortuosity parameter 
of the medium. 

There are obvious disadvantages to this whole procedure in 
defining a Reynolds number. First, a length parameter involving the 
product of the permeability and some vague tortuosity parameter cannot 
be conceived as ‘a physically identifiable length. Second, the medium 
must undergo a series of experiments before this pseudo-length parameter, 
Kb, can be determined. Third, even the most subtle difference between 
the internal geometries of two cores can make a unified correlation 
between them based on a Reynolds number defined by (25) impossible. For 
example, Beavers and Sparrow found that the presence of loose fiber ends 
within a medium made impossible any correlation with a medium that did 
not have loose fiber ends. 

Nevertheless , the proper length parameter is the one which cor- 
relates friction and heat transfer data into single curves, and this is 
the ultimate test. . 
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B. FRICTION FACTOR 



The friction factor is the dimensionless variable which correlates 
the pressure loss behavior of similar surfaces. It is loosely defined 
as the ratio of dissipative forces to inertial forces. Any force (and 
not exclusively viscous forces) which causes a pressure drop across a 
medium is a dissipative force. 

Again there is more than one approach,. The first is the solution 
for the familiar Fanning friction factor, f, based on a channel analysis. 
The second is a friction factor which is an outgrowth of the definition 
of the Reynolds number in terms of the permeability, C^. 

1. Fanning Friction Factor, f : 

Consider the following idealized situation [13]: 
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The pressure drop at entrance consists of two parts. First 
there is a pressure drop due to the decrease in the flow area, which 
from continuity results in a flow acceleration. 

Note that the identity equals the porosity (2) is strictly 

true only for the model which is porous to flow in only one direction. 
For the other two approaches this identity slightly overestimates 
since for these approaches the medium is porous to flow in three dimen- 
sions. 

Second, there is a pressure loss caused by boundary layer 
separation, due to the irreversibilities of an abrupt contraction, and 
manifested by a vena contracta. 
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where is the empirical contraction coefficient, 
loss at entrance is therefore: 



The total pressure 
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Similarly there is a pressure change at exit consisting of two 
First, there is a pressure rise due to the increase in area: 
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Second there is a pressure loss due to the irreversibilities of a 
sudden expansion 

z 



p = K. — 



( 30 ) 



where is the empirical expansion coefficient. The total pressure 
change at exit is the sum of (29) and (30): 
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The remainder of the pressure drop across the heat exchanger is 
due to the overall fluid acceleration and the core friction. The pres- 



sure loss due to overall fluid acceleration is given by: 

Q ■ u ex') - cior^s-r. 
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The pressure loss due to core friction defines the friction factor: 
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After replacement of density by specific volume, the total 
pressure drop across the heat exchanger is : 
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For the isothermal flow of a perfect gas, the pressure is inversely 
proportional to the specific volume and (37) may be solved for f in 
terms of pressures: 
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For heat exchangers where there are successive contractions and 
expansions (such as fibrous materials) the entrance and exit coefficients 
are effectively lumped into the core friction factor [13], Therefore 
in this study: 
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For small pressure differentials, the first term in (39) is by 
far the greatest contributor to the friction factor. A further approxi- 
mation then is: 



= pm= P, -P. 
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Equation (40) defines the Fanning friction factor for fluid flow in porous 
fibrous media. 
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2. Correlating Friction Factor > C^ : 



Since the friction factor is the ratio of dissipative to inertial 
forces, the common practice in porous media studies is to divide the 
pressure gradient as given by Reynolds equation by one-half of the 
inertial term and define the quotient as a correlating friction factor 
[3,9]: ^ 

Cr = 






3^ 



(41) 



Recall that the Reynolds number in terms of the permeability and tortu- 
osity parameter is: 



Np, = 



Substitution of Reynolds equation 



dx - ^ 
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into (41) yields 
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The kinematic viscosity is V - . Thus 
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Cr ~ ^ UKb 1 



(43) 
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It should be noted that ~ 



and that therefore 
equation (41) for might be written .(with ^ replaced by ^P/l ) 
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The term in braces is the first term in brackets in the expression for 
the Fanning friction factor, f, equation (40). Therefore the correlating 
friction factor is related to the Fanning factor by: 
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C. HEAT TRANSFER 

The parameters which characterize the ability of a surface to transfer 

heat include the heat transfer power per unit heat transfer area (h), and 

two non-dimensional forms of h: the number of transfer units (N_^ ) and 

tu 

the Colburn j -factor. The heat transfer power per unit heat transfer 
area, h, is actually the experimental constant in Newton’s Law of Cooling: 

= h A (aT ) (47) 

and is more commonly termed the "convective heat transfer coefficient." 

The number of transfer units, N_^ , is the ratio of heat transferred by 
the surface per unit of time to the thermal capacity rate of the fluid: 

"" ^»A/(rriCp) (48) 

The Colburn j-factor is a heat transfer parameter based on Reynolds 
analogy and empirical observation which is found to behave in a nearly 
identical manner as the friction factor, for most fully-developed internal 
flows [6]: 
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Due to instrumentation difficulties inherent with complex geometries, 
the usual steady state method of determining these characteristic param- 
eters are generally unsatisfactory for compact heat exchanger surfaces, 
and it is necessary to resort to transient techniques. In this ”first- 
cut” analysis of heat transfer from fibrous metals, the simplest such 
technique, the "single blow transient test" or "maximum slope" technique 
was used. This method has been amply described in references [2, 5, 11, 
15, 16, 17, 19, 23, 24, 25, 29, 30, and 31]. 

The problem of analytically determining the fluid and solid temper- 
atures as a function of position and time for a porous solid with constant 
initial temperature and step change in the fluid entering temperature was 
studied independently by Anzelius, Nusselt, Hausen, and Schumann between 
1926 and 1929. The Anzelius solution is overly simplified, while 
Schumann’s solution is valid only for liquids. 

The approach of Hausen, however, appears to be valid for the problem 
at hand. If the fluid is flowing under the Influence of a pressure gra- 
dient in one direction only, the problem may be lumped in the other two 
coordinate directions, and a one-dimensional energy balance may be made 
on a differential element of the solid. 
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Assumptions made are : 



a) Properties of the fluid are* temperature-independent. 



b) Fluid flow is steady. 



c) Porous solid is homogeneous. 



d) Thermal conductivities of both fluid and solid are infinite 



perpendicular to the fluid flow direction, 

e) Thermal conductivity of fluid is zero in the flow direction. 
The boundary conditions are: 

a) The matrix is initially at uniform temperature. 

b) At time equals zero there is a step change in the entering 
fluid temperature. 

c) Matrix boundaries are adiabatic. 

The energy terms (on a per unit time basis) are then: 



a) 

b) Heat transferred to solid by convection = ladV^f “ ^5 

c) Heat transferred from the fluid 





d) Heat transferred from solid by conduction 




Therefore energy balances on the solid and fluid are: 




(50) 




(51) 



The following non-dimensional variables are now introduced: 



b) Generalized time. 



a) Reduced length. 
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c) Longitudinal conduction parameter 






CpU 



Since the thermal capacity of the fluid within the solid matrix 

at any instant of time is much less than the thermal capacity of the 
solid, the second term in the definition of is negligible compared to 
the first . Therefore , 



T = 
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(52) 


Z, andT^ equations (50) and 


(51) become; 
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For the case of zero longitudinal conduction (i.e., 0), the above 

equations become simply: 

: = t,f Ls 

(55) 
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— = t -t 

This system of differential equations, subject to the initial and 
boundary conditions previously specified may be solved analytically 
using Bessel functions. Theoretical fluid temperature - time curves 
at constant Z (in particular at the matrix exit, where ^ ^tu ) 

may be plotted. To determine the of a heat exchanger surface, the 
experimental temperature- time history at the matrix exit is matched to 
one of the family of theoretical curves. This procedure is not only 
tedious , but also depends on the assumption of zero longitudinal con- 
duction. 
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In 1950 Locke [19] eliminated the necessity of matching experi- 
mental and theoretical curves. A new time variable is defined: 






= T/n 



■Lu 



( 57 ) 



He differentiated the analytical solution to Hausen's equations (retain- 
ing the assumption X= 0) and evaluated the derivative at the matrix exit 



(Z = ): 

,tu 









(58) 



where t^ = fluid exit temperature for yU > 0, t^ = fluid exit temper- 
2 ^ 

ature for /x ^ 0, and = steady state temperature. Locke showed that 






as long as the change in the fluid entry temperature is a step change 
(and longitudinal conduction is zero), the of a heat exchanger sur- 
face is uniquely related to the maximum slope of the temperature-time 
curve (i.e., the maximum value of (58)) at the matrix exit, hence the 
name "maximum slope technique." 

Although the assumption of zero longitudinal conduction is essen- 
tially valid at large flow rates (i.e., ‘high Reynolds numbers), it is 
quite invalid at low Reynolds numbers for most surfaces. If these 
longitudinal conduction effects are included in the analysis, the 
maximum slope is a function of A as well as N . With finite differ- 
ence techniques, Howard [11] has included these effects and has tabu- 
lated (Table II) and plotted <J'ig5. 1,2) maximum slope as a function of 
,A). However, for fibrous metals such as steel wool, it is felt 
that any possible longitudinal conduction path would be too lengthy 
for there to be any pronounced effect. Thus in this study Howard's 
results are used only for the special case of A = 0. 
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Kohlmayr [15, 17] showed that the maximum slope technique may be 
inaccurate for 3 . 0 > 2. 0 , unstable for and is singular 

j 

at =2.0 (Fig. 3). He proposed an indirect curve matching technique 
based on the first moments of the areas under the experimental and theo- 
retical temperature-time curves to circumvent this difficulty [16]. More- 
over, this "centroid technique" is not restricted to step changes in the 
fluid entry temperature. As will be mentioned in the Experimental Tech- 
niques section, in this study the step change is approximated by securing 
the electrical power to the heaters in the flow stream. The actual 
response of the air is exponential rather than stepped and therefore 
Kohlmayr’ s centroid technique could be used to correct this deviation 
from Locke’s idealizations. However, Bruno [5] has studied the various 
transient testing techniques and reports that there is no measurable 
improvement in the results by using the centroid techniques for 
Due to the expected large surface area and restricted flow rate through 

fibrous metals, it was felt that all values would exceed 3.5, there- 

tu 

by eliminating the need to correct for deviation from step. 

Thus the heat transfer parameters ^ for fibrous metals were 

determined in this study by Locke’s maximum slope technique, subject to 
the assumptions of exact step change in the fluid temperature, and zero 
longitudinal conduction . 



36 



III. DESCRIPTION OF SURFACES 



Steel wool was selected as a representative fibrous metal. Steel 
wool is low carbon bessemer wire with high tensile strength. It may 



be a continuous piece up to 100,000 feet long. It usually has three 
edges, but may have five l4]. The federal government has set sped- 
fications on the fiber size and type of steel for seven grades between 
0000 and 3. The restrictions on fiber size are listed in Table I [7]. 



Table I. Thickness requirements for individual fibers 



Class 


Not more than 6 
percent under 


Mean width 
of fiber 


Not more than 6 
percent over 


No fiber 
to exceed 




Inch 


Inch 


Inch 


Inch 


0000 


— 


0.0006—0.0010 


0.003 


0.006 


000 


— 


.000&— .0013 


.004 


.008 


00 


— 


.0012— .0022 


.005 


.010 


0 


0.0003 


.0016— .0024 


.007 


.012 


1 


0.0005 


.0018— .0026 


.010 


.016 


2 


0.001 


.0040— .0060 


. .014 ' 


.025 


3 


.002 


.0070— .0096 


.018 


.040 



Grades 0000 through. 0 must be low carbon steel. Grades 1 through 3 may 



be low carbon steel or AIST 430 stainless steel. 

Grades 00, 0, 1, and 3 were selected for testing. In the case 
of Grades 1 and 3 the type of steel was determined by a simple test. 
The plain carbon steel will take the copper out of a CuSO^ solution, 
while AlSl 430 stainless steel will not. All grades tested were plain 
carbon steel. Sample fibers of the four grades were examined under a 
microscope with a movable platform. With the cross-hairs in the eye- 



piece and the platform movement control dial (graduated in hundredths 
of a millimeter) the fibers were checked against the federal specifi- 
cations. Cross-sections of each fiber were examined to determine if 
the samples were three or five sided. All samples indeed were three 
sided. In particular, the samples of Grades 00 and 1 (which were 



37 



manufactured by the Brillo Corporation) appeared to have a very nearly 
equilateral cross-section. 

Such an equilateral cross-section was assumed for all grades tested, 
with side length equal to the average of- the range of mean widths of 
the fiber quoted in the federal specifications. (Table I). The volume 
of the material in each sample is equal to the weight of the sample 
divided by the density of the steel. The assumed cross-section and 
volume of steel specify the heat transfer area, A. (See Appendix B). 

The steel wool was tested at essentially the same porosity as that 
with which it is packaged. The porosity may be crudely estimated from 
article 5.1.1 of [7] which states that one pound of steel wool shall be 
packaged in a sleeve . approximately four inches in diameter and fourteen 
to sixteen inches long. Based on a density for plain carbon steel of 

3 

490 Ibm/ft [18], the porosity should be between 0.979 and 0.983. The 
porosity may be more precisely determined from the weight of the sample 
tested and the volume of the test cell in which it was packed. The free 
flow area is approximately determined from the porosity. 

Complete geometrical descriptions of each sample, based on the 
assumptions stated, are found in Figs. 13, 15, 17, 6 19. 

A core known as Solar No. 4 tested by Bruno [5] was also tested to 
insure the reliability of the reconstructed apparatus. 
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IV, EXPERIMENTAL TECHNIQUES 



The single blow transient test apparatus (Figs. 4,5) was designed 
to conform to the assumptions made for the analytical solution of the 
heat transfer problem (equations (50) and (51)). These idealizations 
are: 

1. The velocity and temperature profiles of the fluid 
flow at any cross-section in the matrix are both 
steady and uniform. 

2. The thermal conductivity of the matrix is finite 
(or zero) in the direction of the fluid flow and 
infinite in the direction normal to the flow. 

3. The thermal capacity of the matrix is much greater 
than that of the fluid contained within it. 

4. The thermal properties of the fluid and matrix 
are both constant and uniform. 

5. The convective heat transfer coefficient is some 
suitable average and remains constant. 

6. A step change in the inlet fluid temperature occurs 
at real time equals zero. 

A specially designed nozzle, flow-straightening screens, and an even 
distribution of heater wires across the channel insure that the first 
assumption (i.e., uniform velocity and temperature profiles) is satisfied. 
Plersall [23] made velocity and temperature traverses at the matrix inlet 
section over the entire flow range and concluded that the velocity 
profiles were indeed uniform and that the temperature profile varied 
by no more than +_ 0.5®F. Wheeler [31] has demonstrated the importance 
of the uniformity of these profiles. 

For fibrous metals it is assumed that there is no effective solid 
area, A^ , available for longitudinal conduction. This has the same 
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effect as assuming zero matrix conductivity in the flow direction and 
therefore satisfies the second assumption. 

As long as the fluid is restricted to gases, the assumption of 
negligible fluid thermal capacity compared to matrix thermal capacity 
is justified. It must be noted that the testing procedures and results 
of this study are not applicable to liquids. 

By limiting the temperature rise to 20°F the assumption of constant 
properties is met. The worst offender is the viscosity, which may vary 
by as much as 1.5 percent from its value at the mean temperature. The 
limited temperature rise also insures that the convective coefficient, 
h, is effectively constant. 

As has been discussed in Section II, the assumption of a step change 
in the fluid inlet temperature cannot be precisely satisfied. In the 
current testing method the step change is simulated by turning off the 
heaters. This actually leads to an exponential decay in the air temper- 
ature. Techniques are available to compensate for this deviation from 
step change, but Bruno [5] has indicated that these corrections are not 
required for ”> 3.5. Since it was expected that the range 
experienced would exceed this value, no effort was made to so compensate. 

The data required for the computation of the various heat transfer 
and flow friction parameters is as follows : 

“ atmospheric pressure (inches Hg) 

- orifice static pressure (inches H^O) 

- pressure drop across orifice (inches H^O) 

A Pm - pressure drop across matrix (inches 

P^ - static pressure at matrix inlet (inches H^O) 

t^ - fluid temperature at orifice (millivolts) 
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d - orifice diameter (inches) 
o 

^ - ratio of orifice diameter to pipe diameter 

CS - chart speed (millimeters /second) 
tg“t^ - downstream temperature response (mm) 

Pressures were measured with water manometers and checked with 
inclined gages where the flow rates and orifice plates permitted. An 
exception was the atmospheric pressure, which was measured with a mer- 
cury barometer. 

Only the matrix downstream temperature response, (t^-t^) , was 
measured, although the upstream temperature monitored to 

check for the 20^F temperature rise from the heaters. A more complete 
description of the equipment and instrumentation is furnished in 
Appendix A. A sample output is shown in Fig. 6. 

Because the turbocompressor which induces the draft through the 
apparatus, would not function for much more than seven minutes with- 
out overheating and tripping the circuit breaker, and because of the 
importance of isothermal flow to the permeability studies, the friction 
factor data runs were made separately from the heat transfer data runs. 
This was in contrast to previous studies where the friction factor data 
was collected immediately after the heat transfer data. 

For a heat transfer data run the pressure drop across the ASME 
standard orifice plate corresponding to the desired mass flow rate 
was predetermined. The number of heaters required to achieve a 20^ F 
temperature rise at that flow rate was also predetermined. When the 
desired pressure • drop had been obtained, the recorder pens were zeroed, 
and the heaters were energized. The channel indicated imme- 

diately whether the 20^F rise has occurred, \flien the core had reached 
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a steady state temperature (as indicated by the trace) all 

pressure measurements and the orifice temperature were recorded. The 
heaters were then de-energized and the transient temperature response 
downstream of the matrix recorded. After the completion 

of the run the maximum slope was visually determined with a straight 
edge from the recorder output and entered in the data sheet. The data 
sheet conformed to the input format for a digital computer program 
written by Trost [30l and modified for fibrous media, which reduced 
the data, (See Appendix B). * 

A friction factor data run has made in the same manner as a heat 
transfer data run except that no heaters and therefore no recorder 
were required. 
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V, PRESENTATION OF RESULTS 



For each of the four grades of steel vrool tested, the heat transfer 
and flow friction characteristics have been computed and appear in 
Tables III to X. Although Fanning friction factors were computed 
for the heat transfer data runs (Tables VII to X), the flow state at 
the time of measurement cannot be guaranteed to have been isothermal 
(particularly at low flow rates). For this reason the friction factor 
results in Tables III to VI are preferred, since these runs were made 
without, the heaters and are the basis for the permeability analysis. 

The Colburn j factor is plotted as a function of N^^^ for each sample. 

The values of N and h as a function of N^ appear in Tables VII to X, 

Figure 7 is a plot of the friction data in the form of (23), Figure 
8 is the same plot in the form of (13), From these plots the sample 
permeabilities and tortuosity parameters are determined. This informa- 
tion is prerequisite for the selection of a length parameter for the 
Reynolds number and therefore appears first in the results. Figures 9, 
10, and 11 plot the Fanning friction factor, f, as a function of N^^^ 
based on 'the three different approaches attempted. The correlating 

friction factor, C-, is plotted against N based on Kb in Figure 12. 

t Re 

Figures 13, 15, 17, and 19 contain the geometric and physical char- 
acteristics of the four samples. Figures 14, 16, 18, and 20 graphically 
present the heat transfer and flow friction characteristics (i.e,, j and 
f versus N^^^ based on Kb) for each grade of steel wool. Figure 21 
compares the j values of the various cores as a function of N^^^ based 
on Kb, Figure 22 compares the heat transfer power per unit of area of 
each core as a function of the friction power per unit of area. 
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VI. DISCUSSION OF RESULTS 



A. RELIABILITY OF APPARATUS 

The Solar No. 4 core previously tested by Bruno was tested in the 
reconstructed apparatus to insure its reliability. An unsatisfactory 
amount of electrical noise was present in the thermocouple signals, 
apparently due to a combination of the unshielded 0.001 inch thermo- 
couples (Appendix A), ultra-sensitive Brush amplifier, Variac power 
supply, and Physics Department linear accelerator, located immediately 
below the apparatus. With this system, N^^ values were found to deviate 
by as much as 23 percent and generally 12 percent from Bruno’s plot. 

The installation of a more sensitive reference junction (see Appendix A) 
may be part of the explanation. Friction factor results deviated by 
as much as 12 percent and were consistently 7 to 8 percent higher than 
Bruno’s plot. 

Since the uncertainties were +_ 7 to 9 percent in N^^ values, 2 
percent in N^^^ values, and 3 percent in f values for Bruno [2, 5] and 
somewhat higher in the reconstructed system due to the previously 
mentioned difficulties (see section on experimental uncertainties), it 
was felt that these results constituted satisfactory agreement. In 
other words, the problem was not felt to be with the thermocouples, 
and the research proceeded. 

Subsequently a less sensitive amplifier-recorder system was installed, 
the Variac was replaced by a motor-generator DC source , and testing was 
performed when the linear accelerator was not in operation. Also a more 
elaborate manometer system was installed (Appendix A). 
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B. PERMEABILITY ANALYSIS 



The permeabilities and tortuosity parameters for the four samples 
were determined by means of a first order polynomial least squares curve 
through the data plotted in the forms suggested by both (13) and (23). 
The former equation assumes isothermal incompressible flow with a con- 
stant pressure gradient. This technique was used by Beavers and Sparrow 
in their study of water flow through porous, fibrous media [3]. The 
latter equation assumes Isothermal flow of a perfect gas, but is not 
restricted to incompressible flow or constant pressure gradients. The 
permeabilities computed by the second method are on the order of one 
percent higher than those computed by the first method for the four 
samples examined. Therefore, compressibility effects are considered 
to be negligible in the computation of permeabilities and tortuosity 
parameters. This is not unexpected, since the maximum change in the 
pressure across the cores tested did not exceed 50 inches of water. 

In Fig. 7 the results of the second type of analysis (i.e., com- 
pressibility effects included) are plotted. With the possible exception 
of grade 1, a straight line plot is a reasonable curve through the data 
points. Outside of the data points corresponding to very small flow 
rates (where the relative error in the pressure measurements is large) 
and those corresponding to very large flow rates (where the absolute 
error is large, due to fluctuations in the flow), the maximum deviation 
of the data points from the linear least squares plot is 11 percent. 

This is not nearly as good as the 2.5 percent maximum deviation reported 
by Beavers and Sparrow in [3]. However, their study was for the flow 
of an incompressible fluid (water) over a velocity span of 20:1, with 
14 pressure taps, cathetometer , calming reservoirs, and direct flow 
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measurement. In this study the fluid was compressible, the velocity 
span was nearly 200:1, with indirect flow measurement (orifice plates), 

2 taps, water manometers read by eye, and no calming reservoirs. With 
all this taken into consideration. Fig. 7 might be considered to be in 
adequate agreement with the empirical Reynolds equation (10). 

Conversely, the data points show a marked concavity rather than 
random deviation from the linear least squares plot. This would indi- 
cate the presence of higher order terms than those in (10). Recall 
that Forchheimer proposed the addition of a cubic term to Reynolds 
equation, but that the coefficient was usually so small that the term 
was negligible except at very large velocities. In light of this and 
the extensive velocity range encountered in this study, a second order 
polynomial least squares curve was passed through the data points plot- 
ted according to the manner suggested by equation (13). (Hereafter 
compressibility effects are ignored. ) The results appear in Fig. 8 
along with the linear least squares curve. In general the data points 
deviate from the. second order least squares curve by less than half 
the amount of deviation from the linear curve. This is considered to 
be strong justification for the replacement of Reynolds equation (10) by 
Forchheimer 's equation: 



The second order polynomial least squares curves plotted in Fig. 8 are 
representations of (.59) recast in the following form: 




(59) 




( 60 ) 
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where compressibility effects have been neglected and 



( 61 ) 

(62) 

(63) 



* 

C = 





(64) 

(65) 



It is not clear what physical significance this additional term has, 

but it evidently involves a coupling of inertial and viscous effects 

and is probably dependent on the fluid. Dimensional analysis shows 

3 6 

that c- has units of (force-time /length ) and that c is dimensionless. 



The particular results from Fig. 8 are: 



Grade 


Permeability 


Tortuosity 




Steel 


K 


Parameter, b 




Wool 


(ft ) 


(i/ft) 


c 


00 


0.479 X IqT'I 


416.8 


-9.57 X lO' 


0 


0.718 X 10 ' 


243.1 


-3.02 X 10' 


1 


1.50 X 10~^ 


201.4 


-3.98 X 10 


3 


6.51 X 10“' 


61.51 


-0.512X 10 



From these results the characteristic length parameters (Kb) are 
0.0000199 ft., 0.0000167 ft., 0.0000302 ft., and 0.0000346 ft. for 
Grades 00, 0, 1, and 3 respectively. 
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The presence of this cubic term raises questions regarding the 
derivation of the Reynolds number based on the product of the permea- 
bility and a tortuosity parameter (25) where one term was said to embody 
all the viscous effects while the other embodied all the inertial effects. 
It is noK evident that there is a third term which may involve both 
effects. Nevertheless, the whole procedure was an approximate one, 
and if the resulting Reynolds number correlates the data there is no 
point in questioning the validity of the derivation. 

C. CORRELATING LENGTH PARAMETER AND FRICTION FACTOR 

The Fanning friction factor, f, was chosen as the test for the 
proper length parameter to be used in Reynolds number. Once the hydro- 
dynamic correlation was obtained, it was thought that the heat transfer 
correlation would involve the same length parameter. In Fig. 9 a cor- 
relation of f and N^^^ is attempted based on a flow passage hydraulic 
diameter from the model porous to flow in one direction. In Fig. 10 
the same correlation is attempted based on the crossed-rod matrices 
model. Figure 11 is the same correlation based on the product Kb. 

It is evident that the second model is better than the first in 

that the curves are somewhat closer together. However the plot of the 

/ 

curves f versus N^^ based on Kb represents the tightest grouping. In 

particular note that the curves for Grades 00 and 0 practically coincide. 

These two grades had practically the same characteristic length (1.99 x 
-5 -5 

10 ft. and 1.67 X 10 ft. respectively). The same is true of the 

curves for Grades 1 and 3 which also have comparable characteristic 

-5 -5 

lengths (3.02 x 10 ft. and 3.46 x 10 ft. respectively). This is 

considered to be a clear indication of the superior correlation to be 

had by using a characteristic length equal to the product of the medium 
permeability, K, and the medium tortuosity parameter, b. 
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The fact that a perfect correlation was not obtained in Fig. ’ll 
may be partly explained by the inverse proportionality between f and 
the heat transfer area, A. This area has the largest uncertainty of 
any variable in the study (see section on experimental uncertainties). 

The correlating friction factor, C^, is plotted against in Fig. 

12. As was demonstrated in the section on theory, if Reynolds equation 
ClO) is valid then is automatically determined by the Reynolds number 
based on Kb: 

The data points show adequate agreement with (44) for Reynolds numbers 
less than one. However, for Reynolds numbers greater than unity, equa- 
tion (44) is seen to over-estimate as defined by equation (41). This 
is no surprise, since there appears to be a cubic term in the empirical 
equation for the pressure gradient. This ciibic term is insignificant 
except at large velocities, which explains the deviation of the data 
points from the theoretical curve in Fig. *12 at high Reynolds numbers. 

D. HEAT TRANSFER AND FLOW FRICTION CHARACTERISTICS 

The heat transfer (j) and flow friction (f) characteristics of 
Grades 00, 0, 1, and 3 are presented in Figs. 14, 16, 18, and 20 respec- 
tively and in Tables III to X. Table VII shows that for Grade 00 when 
j <^0.00675, the corresponding value of N^^ is less than 3.0, where the 
maximum slope method is highly unreliable. Similarly for Grade 0 when 
j<io.0118 and Grade 3 when j< 0.0579. All heat transfer results which 
fall in these ranges must be disregarded, due to the high uncertainty 
in the N_ values. 
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Also the. (t -t ) traces for all but the last five sets of Grade 00 

O J. 

data had an excessive amount of noise, and the determination of the 
maximum slope was at best a crude process. Inasmuch as the "noise” 
died off as the floK rate increased, it is thought that it is actually 
the random fluctuations in the temperature, similar to the fluctuations 
observed when the test section vras removed and only the quiescent ambient 
air temperature was recorded. This is perhaps indicative of some flow 
disorder downstream of the Grade 00 sample, but there is insufficient 
information to draw any firm conclusion. Thus all but the data points 
taken at the five largest flow rates for Grade 00 must be regarded as 
qualitative only. 

The Colburn j-factor is the result of an empirical correction to 
the Reynolds analogy for non-unit Prandtl numbers [6]. The Reynolds 
analogy assumes that for fully developed flow (particularly turbulent 
internal flows), the hydrodynamic and thermal boundary layers are 
similar. Therefore the heat flux at a surface is proportional to the 
shear stress at that same surface for a given fluid [8]. • Thus the 
Stanton number is proportional to the friction factor and the curves 
of f and j as functions of Reynolds number are parallel. 

In previous studies of compact heat exchanger surfaces , the flow 
passages were well defined channels, and the parallel behavior of the 
f and j factors predicted by Reynolds analogy (with Colburn’s correction) 
was clearly observed [29]. This was not found to be the case for any 
of the four samples tested (.Figs. 14, 16, 18, 6 20). Whereas the 
friction factor was found to decrease with increasing Reynolds number, 
the j factor was found to rise to a maximum and then taper off. Figure 
21 indicates that -this peak in the j factor occurs between Reynolds 
numbers (based on Kb) of 0.5 and 1.5 
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There are three possible explanations for this anomalous behavior: 

1. Experimental error 

2. Longitudinal conduction effects 

3* Thermal boundary layer does not behave in the 

same fashion as the hydrodynamic boundary layer. 

If the cause were experimental error, then it would have to be a 

systematic error, as all the cores behave in a similar manner. Trost 

[30] observed the same sort of behavior for Solar No. 1, Solar No. 6, 

Stainless Steel Plate-Fin reference matrix, and the entire series of 

brass cores (with varying L/r^ ratio). He could offer no explanation, 

n 

and none can be offered here from a systematic error viewpoint. It 
might also be noted that no such behavior (i.e., a peaking of the j 
curve) was observed during the testing of the Solar No. 4 core. This 
would seem to rule out systematic error as a cause. 

The inclusion of non-zero longitudinal conduction would tend to 
shift the low flow rate portion of the j curve up to where it would 
be more nearly parallel to the f curve, as demonstrated by Ball [2]. 

It has already been established that an improper assumption of X = 0 
will cause siibstantial error in the heat transfer results at low flow 
rates. However, the small cross-section to surface area ratio and the 
sinuous path of the fibers make it difficult to conceive of such a 
mechanism at work in fibrous materials. A simple experiment to test 
this conclusion is discussed in the section on suggestions for further 
investigation . 

The third explanation is the most plausible. When the thermal 
boundary layer no' longer behaves analagously to the hydrodynamic 
boundary layer, the Reynolds analogy is clearly invalid. This fre- 
quently happens in -external flows, where the turbulence induced by 
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the body contributes to the drag on the body, but only the skin 
friction is important for the transfer of heat. Significantly, 

Colburn [6] noted that "the use of a friction line for predicting 
heat transfer would be very unsafe for flow across tubes." In par-- 
ticular, when separation occurs, it is well known that the j and f 
factors no longer behave in the same manner [18]. Inasmuch as effects 
similar to separation have been postulated as the cause of non-Darcy 
flow [3, 9, 12, 26, 6 27] and inasmuch as virtually all the data runs 
were in the non-Darcy regime, the occurrence of separation or similar 
phenomenom would seem to be the most probable explanation for the 
failure of the f and j curves to behave in a parallel fashion. 

Figure 21 was an -attempt to see how well the j factor for the 
four samples was correlated with a Reynolds number based on Kb. There 
is really too much scatter in the data to draw any firm conclusions. 

The j factor, like the friction factor, is inversely proportional to 
the heat transfer area, A. By the same reasoning as before, if A 
were more precisely known, the curves might be more closely grouped. 

It is of interest to note that Trost [30] experienced the same inability 
to correlate j v/ith (based on a well-defined channel hydraulic 
radius) for six different brass cores. 

Figure 21 also shows that the optimum Reynolds number for steel 
wool as a heat exchanger (neglecting friction considerations) is in 
the vicinity of = 1, since the j factor is a maximum. Therefore 
more heat is transferred per pound of air flow in that range of Reynolds 
numbers, hence the term optimum for heat exchange. 

A figure of merit in evaluating heat exchanger surfaces is the heat 
transfer power at a given friction power level (all evaluated at standard 
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conditions of 500®F and 1 atmosphere), as in Fig. 22. The higher a 
curve is on such a plot, the better the surface it represents is as 
a heat exchanger. Under the conditions of an assumed average equi- 
lateral cross-section for each fiber (from which the total heat 
transfer area is computed), and the assumption of operation at stan- 
dard conditions, Fig. 20 indicates that Grade 3 steel wool is the best 
surface of those tested. It appears that the higher the grade, the 
better it is as a heat exchanger. All grades of steel wool tested 
were evidently better heat exchangers Con a per unit area basis) than 
was Solar No. 4, which was also plotted. 
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VII. EXPERIMENTAL UNCERTAINTIES 



No attempt was made to evaluate errors arising from failure to 
satis:^ the idealizations discussed in Section IV. 

The remaining sources of error are uncertainty in physical con- 
stants, inaccuracy in geometrical measurements, instrumentation error, 
and error in the permeability analysis. The method of Kline and 
McClintock [14l was used to determine the uncertainty in the final 
results. 



A. UNCERTAINTY IN PHYSICAL CONSTANTS 

Values for the physical constants of plain carbon steel were 
obtained from [18], while those for air were from [10]. The uncer- 
tainties (as defined by [14] and different from percentage error) in 
these values appear to be [23] : 

: + 0.5% 






c : + 0.5% 



c : + 0.5% 
P - 

N„ : + 2.0% 

Pr — 



s — 



^ : + 1 . 



5% 



The thermal conductivity, k^, appears only in the longitudinal conduc- 
tion parameter ,"X , which in this study is assumed equal to zero. Since 
the material of construction is specified only as plain carbon steel, 
there is a fairly wide latitude in the value of and in particular, 
c . The latter may be a major cause of inaccuracy [19]. 

o 
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B, GEOMETRICAL MEASUREMENT INACCURACIES 



Contrary to previous studies, this area is a major source of error 
due to the deformable nature of fibrous media and the randomness of the 
internal geometry. The porosity of the sample seems to have the least 
uncertainty, since it vras determined by two independent (though approx- 
imate) methods (see Section III) which agree to within one percent. There 
is variation in the length of the matrix (L) across its cross-section. 

This length also varies with flow rate because a certain amount of de- 
formation occurs as the vacuum downstream grows larger. The free flow 
area, A^, depends on the approach taken to the Reynolds number and the 
uncertainty in it is therefore a function of the accuracy of the approach. 
The crossed-rod matrices model yields a value of A^ approximately 3% less 
than that computed by the other two methods. The heat transfer area, A, 
has the largest uncertainty due to its extreme sensitivity to the assumed 
fiber cross-section. There is an uncertainty of 25% just in the selec- 
tion of a mean fiber width. The assumption of 30% uncertainty in A may 
even be optimistic. However, as noted by London [20] errors in such 
geometrical measurements do not affect the numerical accuracy of the N^^ 
results. The Colburn j factor and convective heat transfer coefficient, 
h, are affected, but the same results may be obtained if future work on 
geometrically similar surfaces is based on the same geometrical assump- 
tions. 

The uncertainties assumed for the various geometrical measurements 



are : 




1 . 0 % 



L : + 3.0% 



A^ : + 1.0% 

fr — 
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A : + 5.0% 
c — 

A : ^ 30.0% 

W : negligible 

s 

(W is the vreight of the matrix), 
s 

C. INSTRUMENTATION ERROR 

Instrumentation errors were primarily due to the pressure measure- 
ments. The uncertainties suggested by Piersall [23] for the various 
pressure measurements were not justified based on experience in this 
study of fibrous media. The following uncertainties are suggested in 
lieu of Piersall *s values: 

P 2.0% 
o — 

AP : + 2.0% 

^P : + 2.5% 

P : + 2.5% 
m — 

P ^ negligible, 
atm ^ ^ 

The orifice temperature was measured with a sheathed copper- 
constantan thermocouple and read by a Leeds and Northrup potentiometer. 
Based on the manufacturer’s calibration of the thermocouple, the esti- 
mated error is one-half of the smallest potentiometer gradation, or 
about 0.0025 millivolt (approximately O.l^^F). The transient tem- 
perature response measured directly from the recorder out- 

put in millimeters. The assumed uncertainty in the maximum slope thus 
measured is + 3.0%, 

Since the N^^. determined is a function solely of the maximum slope 

(for'X= 0), it is appropriate to evaluate the uncertainty at this point. 

The bulk of the N_ values were in the range 22 "^N^ "^5. From Fig. 3 
tu- ^ tu 
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for X = 0 and an uncertainty in the maximum slope of 3%, the uncertainty 

in the N is between 3.3% (N = 22) and 3.75% (N =5). A few values 
tu tu tu 

of were in the vicinity of 2.0, where the uncertainty in the is 
indeterminate. These results must therefore be regarded with skepticism. 

D. PERMEABILITY ANALYSIS ERROR 

Although primary data was used in the determination of the permea- 
bility, K, tortuosity parameter, b, and cubic coefficient, c, of each 
grade, the uncertainty limits cannot be found by the method of Kline 
and McClintock, since these quantities are not the result of single or 
multiple sample experiments, but rather a least squares curve fitting. 

The uncertainties in K, b, and c are a function of the spread (i.e., 
uniform distribution over the mass flow range) and number of data points, 
as well as the accuracy of the individual points. For these reasons 
the relative uncertainties in K, b, and c at 20:1 odds must be estimated: 

K: + 10% 
b: + 10% 
c: 25%. 

E. RESULTS 

The basic formula for the uncertainty w^ in result R(V^, 
at specified odds is: 




where w^, w^, , w^ are the uncertainty intervals at the same specified 

odds for the n variables of which R is a function. Note that w^ has the 
units of R. If the n variables appear in some form in the formula for 
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R as factors, rather than addends, simplification in (66) results 
from dividing through by R. 



Consider, 
mass flow rate 
equation) : 



as an example, the evaluation of the uncertainty in the 
, m. (See Appendix B for the origin of the following 



m 



= (constant) 






'/2 



(67) 



Let P- = P ^ - P . With (66): 
o atm o 



L To V ^ ; 





( 68 ) 



Divide through (68) by (67): 




w* 

The following relative uncertainties (i,e., ^ ) have already 

been postulated (based on 20:1 odds): 

. 0 . 02 , 

'^APo / 



^ 0-00 isr (negligible) 

Therefore the relative uncertainty in the mass flow rate is: 




In a similar manner, the relative error in the mass flux G is found to 
be 1.7%. 
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There is no reason to compute the uncertainty in the Reynolds number 
based on either the uni-dimensional porosity or crossed-rod matrices 
models, as these are completely artificial concepts. Conversely, the 
Reynolds number based on the product Kb does have physical significance 
(i.e. , there is a correct value). The uncertainty for this Reynolds 
number is: 

: • + 14.3%. 

Re — 

The uncertainties in the other results are as follows : 

C^: + 11.4% 
f: + 30.4% 
j : Hh 30.6%. 

It should be noted once again that these uncertainties are based on 
20:1 odds and are not maximum percentages of error. No attempt has been 
made to include the effect of such theoretical assumptions as zero 
longitudinal conduction, step change in temperature, etc., in the 
uncertainty analysis. 
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VIII. CONCLUSIONS 



Bas.ed on the experimental results of this study, the follovfing 
conclusions may be made; 

1. The pressure gradient in fibrous media may be represented by 
a ciibic equation with the bulk velocity as the independent 
variable : 




K is the permeability of the medium, b is a tortuosity 
parameter, and c is a small, negative, dimensionless number 
whose physical significance is unknown. The above equation 
assumes negligible compressibility effects, as found in 
this study. 

2. The product o.f the permeability and the tortuosity parameter 
is superior to the artificial hydraulic diameter determined 
from either the uni- dimensional porosity or cross ed-rod 
matrices models for use as the characteristic length in 

the Reynolds number for flow in fibrous media. The product 
of the permeability and tortuosity parameter is a function 
only of the medium, but is not a conceivable dimension as 
is a hydraulic diameter. Moreover, this quasi-length parameter 
must be determined experimentally, and is very sensitive to 
the subleties in the structure of the fibrous medium. In 
light of the above conclusion that there is a cubic term in 
the equation for the pressure gradient, it is possible that 
the precise characteristic length may involve the dimension- 
less number _c as well. 

3. A friction factor correlation, as a function of N_ based on 
Kb, which is valid for steel wools at N^^^ less than unity is: 

Cf = 

For Nj^^ in excess of unity, this relation over-estimates C^. 

4. The Colburn j factor shows a peak value in the vicinity of 
Reynolds numbers of one for all grades of steel wool. This 
is the optimum operating point in terms of heat transferred 
per pound of air flow. 

5. The thermal and hydrodynamic boundary layers in fibrous media 
do not behave in analagous manners . Separation , which causes 
deviation from the Darcy flovr regime is the probable explana- 
tion for this dissimilar behavior. 
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6. Grade 3 steel wool (by far the largest fiber size) is the 

best’ overall heat exchanger, based on evaluation at standard 
conditions and an assumed equilateral cross-section. In 
general, the larger the fiber size, the better the perfor- 
mance as a heat exchanger. 
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IX, RECOMMENDATIONS FOR FURTHER STUDY 



This first-cut analysis of heat transfer and flow friction charac 
teristics of fibrous media has raised questions which merit further 
investigation. 

1. ' Since the surface area. A, was the largest source of 

uncertainty in the experimental results, the use of gas 
diffusion and other more sophisticated techniques for 
determining specific surface [27] should be pursued. 

2. The possibility of longitudinal conduction should be 
investigated. Since Grades 1, 2, 3 may be either stain- 
less steel or plain carbon steel, and since these two 
metals differ greatly only in thermal conductivities, 

a simple experiment to determine if any longitudinal 
conduction exists is easily performed. 

3. Inasmuch as values of N^^ less than 3.0 were experienced, 
the use of the cyclic testing technique for future studies 
is recommended. 

4. In addition to sampling various grades of steel wool, it 
is suggested that the same grade of steel wool be tested 
at various compactnesses. Thus porosity, as well as fiber 
size becomes an independent variable. 

5. Some further thought must be given to packaging a steel 
wool sample before testing, so that better control over 
homogeneity, dimensions, and boundary conditions may be 
had. 

Prior to any further use of the single-blow transient testing 
apparatus, some improvements in the instrumentation should be made, 

1. Radiation shielding should be applied to all thermocouple 
and heater leads to eliminate noise problems. 

2. Water manometers should be replaced with pressure transducers 
In this way a graphical output is obtained, eliminating the 
need to make a value judgement under the press of time (i.e., 
during the course of a test). 

3. Additional pressure transducers should be Installed down- 
stream of both the matrix and the orifice. These would 
serve as a check, since the inlet pressure plus the 
pressure drop should equal the outlet pressure. 
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4. If all of the above instrumentation improvements are made, 
six pressures would be measured. Since an eight channel 
recorder is available, the orifice temperature could be 
placed on the seventh channel. Thus a complete graphical 
display for each data run would be available. 
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NUMBER OF TRANSFER UNITS VERSUS MAXIMUM SLOPE 
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Figure 1. Versus Maximum Slope with Conduction Parameter 

(Max. Slope, 0.3-2.05). 
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Figure 3. Error in N Relative to an Error in Maximum Slope vs. 
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Test Section 
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Figure 5 . Schematic Diagram of Test Apparatus 
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Figure 6. Sample of Temperature Responses 
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Figure 7. Determination of Permeability and Tortuosity 
Parameter from the Integrated Foma of the 
Polubarinova-Kochina liquation. 
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Figure 8, Determination of Permeability, Tortuosity Parameter, 
and Cubic Coefficient from the Differential Form of 
Reynolds Equation 
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Figure 9 . Fanning Friction Factor as a Function of Reynolds Number 
Based on a Model which is Porous to Flow in One Dimension 
Only. 
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Figure 10. Fanning Friction Factor as a Function of Reynolds Number 
Based on the Crossed-Rod Matrices Model. 
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Figure 11. Fanning Friction Factor as a Function of Reynolds Number 
Based on the Product of the Medium Permeability and 
Tortuosity Parameter (Kb). 
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Figure 12. Correlating Friction Factor as a Function of Reynolds 
Number Based on Kb. 
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Figure 13. Geometrical and Physical Properties of 
Grade 00 Steel Wool 
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Figure 1U. Heat Transfer and Flovj Friction Characteristics of 

Grade 00 Steel ’-.'ool 
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Grade 0 Steel ./ool 
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Figure I 5 . Geometrical 


and Physical 


Properties of 
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0 Steel '.'ool 
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Figxxre 16, Heat Transfer and Flow Friction Characteristics of 

Grade 0 Steel VJool 
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Grade 1 Steel S'ool 
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Figure 17. Geometrical amd Physical Properties of 

Grade 1 Steel 'jool 
80 







Figure 18. Heat Transfer and Flow Friction Characteristics of 

Grade 1 Steel VJool 
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(Divisions in ituti. ) 
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Figure 19: Geomotrical and Physical Properties of 

Grade 3 Steel vJool 



82 



1 




V . 









Figure 20, Heat Transfer and Flow Friction Characteristics of 

Grade 3 Steel Wool 
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Figiire 21, Colburn j -Factor Versus Reynolds Number 
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Figure 22, Heat Transfer Power Versus Flow Friction Power Per Unit Area. 
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GRADE CO STEEL WOOL - FRICTION FACTOR DATA 



MOOT 


NR FR 


FAN F 


CF 


CFTH DIFCF 


12.92 


0.C8 


6.47330 


21.92537 


26.22597 -4.300L0 


31.19 


0.20 


3.29417 


11. 15779 


12.01401 -0.85622 


36. 13 


0.23 


2.89320 


9.79967 


10.64235 -0.84268 


79.19 


0.51 


1.71273 


5.80133 


5.93749 -0.13616 


95.55 


0.61 


1.49829 


5. 07 50 5 


5.26636 -0. 19131 


197.35 


1 .26 


0.93914 


3.18133 


3.58297 -0.40159 


229.65 


1.47 


0. 86304 


2.92372 


3.36308 -0.43936 


262.59 


1.68 


0.80477 


2.72645 


3.19285 -0.46641 


323.11 


2.06 


0.66830 


2.26607 


2.97085 -0.70478 


387.64 


2.47 


0.62147 


2. 10604 


2.30944 -0.70339 


412.33 


2.63 


0.60511 


2.05076 


2.76182 -0.71106 


393. 97 


2.51 


0.64220 


2.17635 


2.79737 -0.62102 


476.94 


3.04 


0.50001 


1. 69489 


2.65802 -0.96313 


5 45 . 47 


3 .48 


0.48459 


1.64306 


2.57531 -0.93225 


623.21 


3.97 


0.46378 


1.57305 


2.50325 -0.93020 



Table III 



GRADE 0 STEEL WOOL - FRICTION FACTOR RUNS 



MOOT 


NR FR 


FAN F 


CF 


CFTH DIFCF 


36. 64 


0.21 


3.16962 


10.52867 


11.71398 -1.13531 


36.55 


0.21 


3.18507 


10.57999 


11.73739 -1. 15740 


50.25 


0.28 


2.60520 


8.65387 


9.08098 -0.42711 


50. 25 


0.28 


2. 57459 


3. 55216 


9.08094 -0.52877 


62.04 


0.35 


2.25236 


7.48183 


7.73305 -0.25122 


62.04 


0.35 


2.24210 


7.44775 


7.73291 -0.28516 


66.88 


0.38 


2. 16408 


7.18859 


7.31540 -0.12681 


66.41 


0.37 


2.17709 


7.23182 


7.35274 -0. 12092 


79. 97 


0.45 


1.81715 


6.03620 


6.44154 -0.40535 


90.90 


0.51 


1 .67776 


5.57323 


5.90757 -0.33435 


108.87 


0.61 


1 .47654 


4.90487 


5.26146 -0.35659 


140. 96 


0.79 


1.25341 


4.16376 


4.51820 -0.35444 


139.21 


0.78 


1.28519 


4.26931 


4.54991 -0.28060 


171.52 


0.97 


1.19471 


3.96889 


4.06849 -0.09961 


218.10 


1.23 


1.03267 


3. 43075 


3.62558 -0.19434 


244.39 


1.38 


0.98195 


3.26238 


3.45068 -0. 18829 


161.00 


0.91 


1.25312 


4.16286 


4.20287 -0.04001 


245.19 


1.39 


0. 97594 


3.24241 


3.44269 -0.20028 


294.26 


1.66 


0.852‘^4 


2.83233 


3.20294 -0.37061 


345.50 . 


1.95 


0.79664 


2.64723 


3.02507 -0.37784 


324.35 


1.83 


0.75581 


2.51144 


3.09513 -0.58369 


414.06 


2.33 


0.69448 


2.30816 


2.85859 -0. 55043 


502. 38 


2.82 


0.68794 


2.28707 


2.70936 -0.42228 


616.67 


3.46 


0.65043 


2. 16335 


2.57736 -0.41401 


457. 57 


2 .57 


0.67642 


2.24841 


2.77816 -0.52976 


585.91 


3.30 


0.59726 


1. 98620 


2.60649 -0.62029 


769.55 


4.32 


0.53332 


1.77524 


2.46311 -0.68787 


676.39 


3 .79 


0.56239 


1.87105 


2.52706 -0.65601 



Table I? 
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I 






RUN 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 



GRADE 


1 STEEL WOOL - 


FRICTION 


FACTOR DATA 


MOOT 


NR FR 


FAN F 


CF 


CFTH 


DIFCF 


10.80 


0.10 


5.70884 


15.97843 


21.17030 


-5. 19187 


12.60 


0.12 


5.94058 


16.62704 


18.42164 


-1.79459 


13.09 


0.13 


5.83343 


16.32715 


17.78107 


-1. 45391 


15.07 


0.15 


5.38046 


15.05933 


15.70058 


-0.64125 


30.26 


0.29 


3.09263 


8.65595 


8.82646 


-0.17050 


34.24 


0.33 


2.89335 


8.09819 


8.02800 


0.07019 


33.78 


0.33 


2.92283 


8.18072 


8.11408 


0.06664 


72.47 


0.70 


1.71666 


4. 80481 


4. 84508 


-0.04027 


70.21 


0.68 


1.78999 


5.01006 


4.94356 


0. 06650 


84. 10 


0.81 


1.49030 


4.17129 


4.45738 


-0.28610 


96.02 


0.93 


1.49104 


4.17339 


4.15975 


0.01364 


161.22 


1 .55 


1 .09987 


3.07866 


3.28708 


-0.20842 


200. 85 


1.94 


0.99511 


2.78555 


3.03164 


-0.24609 


266.17 


2.58 


0.82385 


2. 30638 


2.77595 


-0. 46956 


314.91 


3 .05 


0.77040 


2.15694 


2.65524 


-0.49830 


385.99 


3.75 


0.71334 


1.99751 


2.53396 


-0.53645 


389.46 


3.78 


0.72909 


2.04164 


2.52927 


-0.48763 


430. 19 


4.18 


0.65172 


1.32518 


2 .47866 


-0.65348 


522.58 


5.07 


0. 57303 


1.60526 


2.3943 8 


-0.78911 


525.58 


5.10 


0.56461 


1.53168 


2.39202 


-0.81033 


593.66 


5.77 


0.57760 


1.61851 


2.34665 


-0.72815 


588.20 


5.71 


0.57416 


1.60883 


2.35017 


-0. 74134 


708.37 


6 .89 


0.47077 


1.31983 


2. 29045 


-0. 97062 


804. 55 


7.82 


0.43115 


1.20936 


2.25574 


-1.04638 
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38 

39 
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GRADE 3 STEEL WOOL - FRICTION FACTOR DATA 



MOOT 


NR FR 


FAN F 


CF 


CFTH 


DIFCF 


954.46 


10.57 


0.60511 


1.65333 


2.18923 


-0.53590 


954.34 


10.57 


0.60542 


1.65417 


2.18925 


-C. 53508 


1005.70 


11.15 


0.59015 


1.61302 


2.17945 


-0.56642 


796. 12 


8.83 


0.62457 


1.70486 


2.22645 


-0.52159 


795.26 


8.82 


0.62590 


1.70849 


2.22669 


-0. 51820 


601 . 17 


6.67 


0.67412 


1 . 83 84 7 


2.30004 


-0.46157 


601.29 


6.67 


0.67390 


1. 33788 


2.29998 


-0.46210 


507.78 


5.6 4 


0.73439 


2.00216 


2.35465 


-0. 35249 


508. 42 


5.65 


0.73120 


1.99347 


2.35417 


-0.36071 


424.40 


4.71 


0.74522 


2. 03120 


2.42421 


-0. 39301 


420.49 


4.67 


0.76954 


2.09746 


2.42323 


-0.33078 


330.36 


3.67 


0.79832 


2. 17543 


2.54521 


-0.36978 


325.40 


3.61 


0.83353 


2. 27136 


2.55352 


-0.28216 


261.52 


2 .90 


0.94764 


2.58202 


2.68863 


-0. 10661 


261.53 


2.90 


0. 94753 


2. 58170 


2.63361 


-0.10691 


261.13 


2.90 


0.93404 


2.54496 


2.6 8966 


-0. 14470 


221.73 


2.46 


0.99255 


2.70422 


2.81169 


-0.10747 


221.73 


2.46 


0.99251 


2.70409 


2. 81168 


-0. 10759 


221.72 


2.46 


0.98500 


2.68365 


2.81171 


-0. 12806 


166. 11 


1.85 


1.08964 


2.96907 


3.08341 


-0.11434 


166.16 


1.85 


1.08993 


2.96934 


3.08244 


-0. 11311 


162 . 92 


1 -.81 


1.13378 


3.08878 


3.10397 


-0.01519 


169.83 


1.89 


1.10759 


3.01746 


3.05967 


-0.04221 


169.83 


1.89 


1.10762 


3.01753 


3.05968 


-0. 04215 


167.88 


1.87 


1.18670 


3.23296 


3 .07202 


0. 16094 


167.41 


1.85 


1.06656 


2.90567 


3.07917 


-0.17350 


167.42 


1.85 


1.05316 


2.86917 


3.07914 


-0. 20997 


154.68 


1.71 


1.12328 


3.06016 


3.16907 


-0.10891 


124.85 


1.38 


1. 19706 


3.26107 


3.44936 


-0. 18829 


100.70 


1.11 


1.46836 


4.00009 


3.80077 


0. 19932 


104. 63 


1.16 


1.53966 


4.19431 


3.72464 


0.46967 


68.99 


0.77 


1.73339 


4. 72198 


4. 61429 


0. 10769 


57.90 


0.64 


2.34971 


6.40092 


5.11421 


1. 28670 


41.96 


0.47 


2.13194 


5.80765 


6.29563 


-0.48798 


36. 82 


0.41 


1.93862 


5.28101 


6.89474 


-1. 61373 


29.72 


0.33 


2.12317 


5.78375 


8.07186 


-2.28811 


23.62 


0.26 


2.35244 


6. 40830 


9.63802 


-3.22972 


20.41 


0.23 


2.69709 


7.34714 


10.85471 


-3. 50757 


12.61 


0.14 


4.70579 


12.81903 


16 .35512 


-3.53609 


9.74 


0.11 


5.89741 


16. 06513 


20.59951 


-4. 53439 
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GRADE 00 STEEL WOOL - HEAT TRANSFER DATA 



ooh-msOvoONro‘'(^j:'-vOC7'oo 

^ •#•••••••#•••• 
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00000 000 00 0000 



00 000 00 00 000 00 
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Q ooooooor-Hr-»'Nr\jc\jmvt 



a: 
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< XCL 
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O h - 00 sO rvj O o CO ro rH > 4 “ 
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Table VII 
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GRADE 0 STEEL WOOL - HEAT TRANSFER DATA 
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< 

UJ 

X 



LU OcnUACACOi-t -0 CM Oca rococo ^M’stCMCAvO OfAlA Ola CT'O^ A- LA CM CM O^sO 00 »-H 

xa. OOoOOOCMfAlALAvOLA»-4cMAjA-h-AJCMCMCMC7'Li'^^NtCM(AuO'Or-C7'OOsOOLA>rlA(A 

<0 <^^(^)^tCAU>^•^-U'^OCOCOfMC^>J'4-LA>j“OOrw,-IC^^OOCM^ALACMlAC^OLA^CMCAsO 

X-J OOCMCACAcA 00 CMAr-OOCMCMlA' 0 >t^CMlAC 0 >i“CAa'fAOC 0 .-<^cAOa'' 0 A 400 
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K 
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z 

o 



c^ 
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a: 
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CA h- vO CA O O' LA CM CM ro n O lA LA nO CM CM O' >t 00 ^ h- O' ^ lA CM r- sO 00 

C^ CA LA CA vt 00 ^ CM »H CO lA C< ) ^ CT' o O CA 0> • O sO CM o LA CT' O 00 O CM st LA O si* O vO 
>0 ^ 00 lA si- 0^ Q> M" 00 CA vt >T LA r\J n nO 00 CM M“ O fA CX) rs- O CT' LA O nO ^ 
lA LA — C h- r- fA O fA LA o vf LA 0> >T CA 00 ^ ^ vO o o lA O O A- O 00 CM 00 0> sO >t 
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APPENDIX A: DESCRIPTION OF EQUIPMENT 



A. AIR SUPPLY AND FLOW MEASUREMENT 

A 30 HP, multi-stage Spencer Turbo-Compressor rated at 550 CFM 
Cfor intake at atmospheric pressure) provided an induced draft through 
the test section. 

The air flovr rate was measured with an AS ME standard orifice section. 
The concentric orifices had throat diameters of 0.308, 0.462, 0.775, 
1.232, 1.540, 1.971, and 2.310 inches. These values correspond to 
throat diameter to pipe diameter ratios (|B) of 0.10, 0.15, 0.25162, 

0.40, 0.50, 0.64, and 0.75 respectively. Pressure taps were located 
one pipe diameter upstream and one-half pipe diameter downstream of 
the orifice plate, as specified by [1]. 

B. HEATER SYSTEM 

The heater system used was essentially that used by Bruno [5], 
although it was completely rewired. Fourteen separate bakelite frames 
were each wound with two parallel-connected heater elements. The heater 
elements were 0.0031 inch nichrome vrire, (Fig. 23) which was selected 
largely due to the comparatively low coefficient of thermal expansion 
of nichrome. Each heater element has a nominal resistance of 800 ohms. 

An individual switch controlled both heaters per bakelite frame (i.e., 

14 such switches) and these switches were in turn controlled by a master 
switch which was used to cut off the power to the heaters to simulate 
the step change in air inlet temperature . Thus the number of heaters 
(as well as the voltage supplied) could be adjusted to provide the 20^F 
temperature rise for any flow rate up to approximately 1,000 pounds per 
hour. 
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C. HEATER POWER SUPPLY 



A General Electric Speed Variator motor-generator, operating on 220 
volts AC at 35 amps and producing 250 volts DC at up to 28 amps, provided 
the electrical power to the heaters. The actual voltage supplied was 
regulated with the rheostat on the field winding of the MG set and at 
the test section with a 100 ohm Ohmite wound resistor (for currents less 
than 4 amperes ) . 

The current and voltage at the source were monitored with a Roller- 
Smith 150 a. ammeter and a Simpson multimeter, respectively. At the 
test section the current and voltage were monitored with a Westinghouse 
10 a. ammeter and 300 v. voltmeter. 

D. MATRIX HOLDER AND TEST SECTION 

Both the holder and test section (Fig. 24) were fabricated from 
plexiglass. The matrix holder is a drawer which slides into the test 
section. By removing the rear, -upstream, and downstream faces of the 
holder, the test section of steel wool may be inserted with minimum 
packing effort, thereby causing negligible inhomogeneity in the porosity. 
The flow channel is 3.1875 inches square and samples up to approximately 
3 inches long can be handled. A network of seven .003 inch steel wires 
was strung across the downstream face of the holder to restrain the 
steel wool from being sucked down the duct. These wires were not thought 
to interfere with the friction or heat transfer performance of the sample, 
and experience indicated that they are probably unnecessary anyhow. A 
similar network was placed in the upstream face to prevent interference 
between strands of steel wool and the upstream thermocouples. This net- 
work is also thought to be unnecessary. Styrofoam insulation lined the 
interior of the matrix holder and also served to hold the steel wool in 



place . 
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The test section, into which the holder slides, contains the matrix 
upstream and downstream thermocouples and pressure taps. In previous 
installations the matrix downstream temperature was measured by thermo- 
couples mounted in the holder, rather than the test section. However, 
the possibility of shorts or breakage from pertruding strands of steel 
wool dictated the movement of the thermocouples downstream by approxi- 
mately 0.5 inch. No significant error in the maximum slope technique 
is thought to have been introduced by this movement of the most impor- 
tant sensor. 

E. INLET CONE AND FLOW STRAIGHTENER 

This section was designed and tested by Piersall [23], and was 
found to provide a uniform velocity profile for the air entering the 
matrix. 

Additionally an alteration was made in the radiation shielding of 
the reference thermocouples , located at exit from the inlet cone , which 
may have improved the uniformity of the velocity profile. Originally 
all the reference thermocouples were contained in one 3/8 inch aluminum 
tube mounted across the center of the channel. In this reconstruction 
these thermocouples were placed in five 1/8 inch nickel-plated tubes , 
evenly spaced across the channel. This is thought to delocalize the 
intensity of the disturbance caused by the presence of the thermocouples 
in the stream. 

F. PRESSURE MEASURING SYSTEM 

One pressure tap was located in both the upstream and downstream 
face of the test section (i.e,, on either side of the matrix). As 
previously mentioned, two taps were also located at the orifice section. 
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Each pressure tap was connected by Poly-flo flexible tubing to a high 

pressure or low pressure manifold and protected by a Hoke quick-closing 

valve. The high and low pressure manifolds were connected to the high 

and low pressure sides of the following instruments: 

2 Merriam Instrument Company, 50 in. manometers 
1 Ellison Draft Gage Company, 0-2 in. inclined manometer 
1 Ellison Draft Gage Company, 0-4 in. inclined manometer. 

These manometers were also protected with quick-closing valves. Thus, 

with the proper valve lineup, each of the four pressures measured per 

data run was simultaneously measured on four manometers until the mass 

flow rate and orifice plate combination required that the more limited 

manometers be isolated. 

Atmospheric pressure was measured by a PRINCO mercury barometer. 

G. TEMPERATURE MEASURING SYSTEM 

Temperature sensors were located at the exit from the inlet cone 
(t^), entrance to the test section (t^)) exit from the test section 
(t^), and 15.5 inches upstream from the orifice section 

Temperature t^ was measured by two different sets of five 30 gage 
iron- cons tan tan thermocouples connected in series. Each thermocouple 
was wrapped in teflon tape to prevent shorting. As previously mentioned, 
the reconstructed apparatus has different housing for these thermo- 
couples than that designed by Traister [29]. The ten thermocouples 
were housed two to each of five nickel-plated brass tubes. The 
two thermocouples in each tube were from the two different sets of t^ 
sensors. The tubes were cut-away to the front while closed to the 
downstream side, where they faced the heaters. This system is thought 
to provide superior sensitivity and radiation shielding to the old system, 
where all ten thermocouples were housed in one aluminum tube. 
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Temperature was measured by a set of five .001 inch diameter 
iron- const ant an thermocouples in series. The output of this thermopile 
was referenced against one of the t^ thermopiles so that the net output 
to the recorder was 

Temperature t^ was measured exactly as t. , so the output to the 

o . 

recorder was Ct -t ). 

sj -L 

The two outputs (t^-t^) and Ct^-t^) were amplified by Hewlett- 
Packard 8803 A DC amplifiers and recorded on a Hewlett-Packard 7702B 
two- channel recorder. 

The orifice temperature, t^, was measured by a 30 gage cooper- 
constant an thermocouple sheathed in stainless steel. It was referenced 
to an ice junction and the output read on a Leeds and Northrup Milli- 
volt Potentiometer. 

H. NOISE CONTROL SYSTEM* 

A substantial amount of noise in the (t -t ) and (t^-t ) signals 
was frequently encountered. To filter this out, Mallory 2.0 microfarad 
capacitors were placed between the negative side of the signal and 
ground. The time constant of the circuit with the addition of these 
capacitors was less than the minimum time constant of the heater- thermo- 
couple system reported by Bruno [5]. 

Additionally a Sprague 1000 micro-farad capacitor was placed in the 
power supply circuit to filter any noise from the motor-generator. This 
had no effect on the transient response of the apparatus, since the 
master switch isolated the power supply from the heaters. 



99 



APPENDIX B: DATA REDUCTION RELATIONSHIPS 



A. GEOMETRICAL RELATIONSHIPS 

Once an average cross-section for a fiber is assumed, the heat 
transfer area of the medium is determined. In this study the average 
cross-section was assumed to be an equilateral triangle with side 
length w, where w is the median mean width of the fiber specified 
for the particular grade steel wool by EVU. Therefore: 



Cross-sectional area 
Perimeter of cross-section = 3w 



A... =W^(V3’/4') 

fiber 



Hydraulic diameter of fiber = W j^l 3 

Volume of steel in sample 



(B-1) 



V =(A_., ) (Total length of fibers) 

steel fiber ^ 



Heat transfer area = A = (Perimeter) (Total length of fibers) 



But the volume of steel is equal to the weight of the sample divided by 



the density of steel. Therefore: 



Steel 



= 4-V3"Ws/^5w) 

Porosity =^ = ^ ~ 



Specific surface = , 



Free flow area = ^^"^f^fr crossed rod matrices) 



(B-2) 

(B-3) 

(B-4) 

(B-5) 



Equations (B-1, B-2, B-3, 6 B-5) are substituted into equations (1) through 
(6) as appropriate to determine the geometrical characteristics for use 
in the various models. 
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B. MASS FLOW RATE 



The mass flow rate was calculated from the ASME Power Test Code El] 
with Murdock’s 121] interpolation scheme. 



C = orifice coefficient of discharge 

^ = ratio of orifice diameter to pipe diameter 

d = orifice diameter in inches 
o 

= orifice plate thermal expansion factor 
Y = fluid thermal expansion factor 

^ = specific weight of fluid flowing = P/RT if a perfect gas 

2 

P = absolute static pressure at orifice (ibf/ft ) 

R = specific gas constant for air = 53.35 ft-lbf/lbm-®F 
T = absolute temperature at orifice (®R) 
aP^ = pressure drop across orifice in inches of water 

Substitution of the identities for K and in (B-6) yields: 




(B-6) 



where 




(B-7) 



From Fig. 40A of [1] one finds: 




where 



k = Cp/c^ = 1.4- O..V') 



X = aR/P 



(any consistent units). 
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Also from [1], Fig. 38: 



F = 1 
a 

P = (P ^ - P /13.6)C70.75)lbf/ft^ 

atm o 

^atm ~ atmospheric pressure in inches of Hg 

P^ = static pressure upstream of orifice plate in inches H^O 



When all these identities are substituted into (B-7) with the units 
adjusted to give m in Clbm/hr), the equation becomes: 



S89.S/C 


dXl-(04l - O^srt' 


VI 


99.02 



1 

> 

^ J 


AR(Rim - ^ /\3.(o') 


%.g)_ 


to +459.7 



The orifice coefficient of discharge is found from Fhirdock's tables 
and interpolation scheme. 



C. MAXIMUM SLOPE 

The following sketch of the general time-temperature history down- 
stream of the core will illustrate the relationships involved. 
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m>- 




A 










The maximum slope of the general time -temperature history is a function 

of N . and X . 
tu 



J (f/Nj 






0 



The non-dimensional variables are: 



— generalized time variable ^ UA©/('^s^s'^ 

N. = hA/mCp 

iw. ‘ 



Therefore 



s 



(B-9) 



Furthermore 

t3 - i.- _ 

t. - t: 

and its derivative is 



i.-i. 

t - 1\ 



* 1 



*^1 ^ ‘ ' 



a. - t;) 



(B-10) 



Combination of (B-9) and (B-10) yields: 

*1 

WsCs 






d (ti" 



tnCp (t. - t;) 



(B-11) 



MAX 
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In terms of the recorder output the various quantities in (B-11) are: 






= v/x 

X/CChart speed) = X/C5-ci© 



dCtg tj^) Y 
= m 



Equation (B-ll) becomes in these terms: 

- t;' 



4^) 



dr/Mj 



Wc c 



s '-s. 



irn c. 



Y 






cs 



(B-12) 



This value of the maximum slope and the known value of X(in this study, 
X= 0) are used to enter Table II or Fig. 1 to determine the correspond- 
ing value of N 



tu 



D. COLBURN J FACTOR 

The Colburn j factor is defined as 

Vs 






4- Np, 



Vs 



(B-13) 



In this case the mass flux, G, is based on rather than although 

at 98% porosity there is not much difference. Substitution of G = A/A^ 
and multiplication of (B-13) by A/A yields: 

hA . 

J r^p A 



But N = (hA/mc ) and therefore.*' 
tu p 






V3 

Pf- 
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E. HEAT TRANSFER POWER AND FRICTION POWER 



The heat transfer power per unit area per degree temperature 
difference is : 






(where 1 is the characteristic length used in N ). 

Ke 



At standard conditions of 500°F and one atmosphere [13] 

c = 0.2477 BTU/lbm-°F 
. P 



= 0.0678 Ibm/hr-ft 
yO = 0.0413 Ibm/ft^ 



= 0.671 
Pr 



(B-14) 



Substitution of these values into (B-14) leads to: 

BTU 






hr-PC-"F 



(B-15) 



The flow friction power per unit of area is [13] 

r ■ / . . 3 \ 2> 

‘Re 






(B-16) 



Evaluation of (B-16) at standard conditions yields: 

\3 



-7 



F = 1 . 11^(0 

STD 




loN. 






(B-17) 



F. OTHER RELATIONSHIPS 

The relationships pertaining to the Reynolds number and the two 
friction factors have been completely derived in the section on 
theoretical background. 
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